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Abstract 

We analyze the effect of correlations in a simple model of small world network by obtaining exact 
analytical expressions for the distribution of shortest paths in the network. We enter correlations 
into a simple model with a distinguished site, by taking the random connections to this site from 
an Ising distribution. Our method shows how the transfer matrix technique can be used in the 
new context of small world networks. 
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I. INTRODUCTION 



Real networks, like social networks, neural networks, power-grids, and documents in 
World Wide Web (WWW) 2I 0, U, [| , can be modeled neither by totally random networks 
nor by regular ones (see |y, [7[ and references therein for review). While locally they are 
clustered as in regular networks, remote sites have often the chance of being connected via 
shortcuts, as in random graphs, hence reducing the average distance between sites in the 
network. 

In a regular networks with N vertices, the average shortest path between two vertices < / > 
and the clustering coefficient C scale respectively as < I >~ iV , and C ~ 1. The clustering 
coefficient C is defined as the average ratio of the number of existing connections between 
neighbors of a vertex to the total possible connections among them. In random networks, 



however we have < I >~ logiV and C ~ 1/N 



The properties of many real networks, are a hybrid of these two extremes, that is in these 
networks one has < / >~ logN, and C ~ 1. These two effects called collectively small 
world effect, are attributed respectively to the presence of shortcuts and the many inter- 
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connections that usually exist between the neighboring nodes of such networks ^ 
In 1998, Watts and Strogatz introduced a simple model of network showing the small 
world behavior, which since then has been investigated as a model of interconnections 
in many different contexts, ranging from epidemi olog y [iJ . Q|, to polymer physics 
[17I lisi and evolution and navigation [3, 21, 122], The original model of Watts and 



Strogatz contained a free parameter p, by varying which one could interpolate between ran- 
dom and irregular networks. Their model, called hereafter the WS model, consists of a 
ring of N sites in which each site is connected to its 2k nearest neighbors, hence making 
a regular network. After this stage, each bond is re-wired with probability p to another 
randomly chosen site. The value of p tunes the amount of randomness introduced into the 
network. Since there is a finite probability of disconnecting the whole network in this way, 
Newman and Watts j3] modified the model by replacing the re- wiring stage by just addi- 
tion of shortcuts between randomly chosen sites on the ring. Since then many more variants 
and generalizations of small world networks and their different characteristics ( e.g. their 
topology, the properties of random walks on them, etc.) have been studied. Of particular 
interest are three classes of studies. The first class , in which the static properties of small 
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G. 1: Two configurations which have different weights in our calculations but equal weights in 



world networks have been investigated 



24. 
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26l . |27[,the second class , where dynamical 



aspects have been studied [28, 29] and the third class, in which evolving networks are consid- 
ered |3oi I3I I in order to generate small world networks with various connectivity distributions 
, including scale free distributions. 

In this letter we want to consider another variant of the small world network, one in which 
correlation of neighboring nodes in making connections to remote sites is taken into ac- 
count (ie;the presence of a shortcut between two sites affects other shortcuts in the neigh- 
borhood) . For example a node need not make a short-cut to a remote site if there is such a 
connection in its neighborhood. In such networks then, correlations play an important role. 
However to perform such a study by exact non-mean field methods requires a simplification 
in the original model. We assume all shortcuts are made via a distinguished site at the center 
of the ring. More than being a simplification, this type of network has practical relevance in 
many situations where a central distinguished site governs all the remote interconnections. 
We note in passing that such central sites accommodating a large number of connections, 
may exist either in the architecture of the original networks or else may appear dynamically 



in evolving networks 
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In this way we assume that contrary to the original model, [lj, the 



two configurations in figE] , both with 5 shortcuts are not equiprobable. 
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FIG. 2: A simple model of small world networks. 

II. THE MODEL AND SMALL WORLD QUANTITIES 

We consider a circular network of N vertices with a distinguished central site 121 The 
links on the ring have unit length. Each shortcut connecting any two sites on the ring is 
also of unit length. We assign a random variable Sj G {0, 1} to each site i of the ring. This 
random variable is 1 or according to whether the site is connected to the center or not,fig|21 
Any configuration of these spin variables corresponds to one and only one configuration of 
connections to the center. For example in fig^if each bond is independently connected to 
the center with probability p, then the probabilities of both configurations are equal and 
proportional to p 5 (l — p) u . In general and in the absence of correlations we will have: 

p = ^ ( P \si+sg+-s N /i\ 

Z [ l-p ) ' {) 

where Z is a normalization constant. To consider correlations we generalize the above 
distribution to an Ising type distribution, namely to: 

1 / N \ 

p {^} = ^(n^c^ +i ). (2) 

For ( = 1 and r = we obtain the original model of The value of ( controls the 
correlations. 

First let us consider the directed model, i.e; the links on the circle are directed, say clockwise. 



4 



FIG. 3: A typical configuration which decreases 

Looking at figEU we consider a typical configuration like the one shown in this figure ,m which 
the nearest shortcuts to sites 1 and j are connected at sites i and k. This configuration 
reduces the distance between sites 1 and j by an amount k — i — 1. Not that the sites 
between i and k may or may not be connected to the center. In any such configuration the 
quantity X it k(l,j) defined as: 

(1 - si) • • • (1 - s i _i)s i s fc (l - s fc+1 ) • • • (1 - Sj), (3) 

takes the value l.The average of this quantity gives the probability of such a configuration. In 
order to find the probability of the shortest path between sites 1 and j to be equal to I, we 
have to sum over all those configurations which give such a shortest path . For I ^ j — 1 the 
above probability is given by: 

i 

p(hr,i) = E < x i,j+i-i-i{hj) >, (4) 

i=l 

where we have used < .. > for averaging over configurations . Normalization determines 

pO-J'J - 1) via: 

p(l,j;j-l) = l-$>(U;0- (5) 

1=1 

The probability that the shortest path between two arbitrary vertices be of length I, is 
obtained from: 

N 

p(l)=l/N £ p(l,j;J). (6) 

j=l+l 



Now the average shortest path between two randomly chosen sites is: 

N-l 

<i>=J2W)- (7) 

1=1 

All the above quantities can be calculated by the transfer matrix method, in which we write 
the unnormalized distribution El as product of matrix elements of a matrix T : 



T 



(8) 



with eigenvalues: 



A± = -[l+rC± v/(l-rC) 2 + 4r]. (9) 



The partition function is Z = A+ + A_ and the number of connections per site is given by: 

T (9 

"■=N¥r inZ - < 10) 

We now consider the continuum limit of the lattice, where the number of vertices goes to 
infinity and the lattice constant goes to zero as so that the periphery of the lattice is 
kept constant at 1 . We then set: j/N — > x,k/N — > s,i/N — > t l/N — > z and 
NXi,k(l, j) — > X(t, s)(x) where the explicit form of the function X(t, s)(x) will be deter- 
mined later. We will then have < x, t, s, z < 1. Here x is the distance along the ring. 
Furthermore we take Np(l,j,l) — > Q(x,z) , therefore : 

Q(x,z) = [ X(t,x + t-z)( x)dt, 
Jo 

Q(z) = [ 1 Q(x;z)dx, (11) 

J z 

where Q(x, z)dz is the probability that two points whose distance along the ring is x have 
a shortest distance between z and z + dz. Then Q(z)dz is the probability that the shortest 
path between any two points be between z and z + dz . So Jq dzQ(z) = 1 and finally : 

< z > = f 1 zQ(z)dz. (12) 



A. The scaling limit 

Intuitively we expect that in the scaling limit, when N — > oo, and r = ^, if we keep £ 
finite, then the number of connections to the center remains finite and in an infinite lattice 



the configurations of these connections become quite sparse and hence correlations can not 
play a role, at least to leading order. Exact calculation also verifies this expectation. Here we 
will consider a different scaling limit where N — > oo, and r = ^ while ( = Co^- This means 
that the tendency of an individual one of whose neighbors has been connected to the center, 
depends also on the total population. This assumption is not far from reality, specially 
in cases where the center approves a limited amount of connections and the applicants 
, competing for connections, are aware of this restriction. It turns out that the model shows 
three distinct behavior according to the value of the parameter r ( . 
For r Co > 1 we will have: 

r 



K = r Co + 



AT(roCo-l) : 



and from we find: 



x -^-W^try (13) 



p=l + 0(i), (14) 



which means that the whole lattice is filled with connections. Also for roCo = 1 we obtain 
p = |, which is also far from small world regime. To be in the small world regime, we should 
keep roCo < 1 , which is the case that we will study in detail. 
In this case we have: 

A + = l + r ° 



Wfl-roCo) 

A - = r ° C ° - N(l-rM- (15) 
Also from (fTU|) we find the total number of connections to be the finite value: 

M :=Np= r ° (16) 
(1 - r (o) 2 

To calculate Aj^l, j) we note that since p — » as 4 , the value of these quantities where 
either or both of % and k take the extreme values 1 or j are suppressed. Using the transfer 
matrix technique, we obtain from ((21, EJ) that for 1 < i < k < j, 

< X i>k (l,j) >= T^T 01 (T k -% 1 T 10 Ti - k -\T N ^+ 1 ) 00 = r(T fc -)n(T^ +1 )oo, (17) 

where T% =< i\T\j > m and (T"%- =< i\T m \j >. Diagonalizing T, using (0,13), and taking 
the continuum limit, we find after some algebra: 

X(s - t)(x) = M 2 e- Mx e M ^- l \ (18) 



where M := 1 _^ Q ^ ■ Inserting this value in and integrating we find: 

Q{x;z) = zM 2 e~ Mz , (19) 

Q{z) = z{\ - z)M 2 e~ Mz + (1 + zM)e~ Mz , (20) 
Turning to (fT2|) it is obtained : 

<z>= h [2+e ' M) -Tp {l - e ' M) - (21) 

As expressed in [1], these relations already hint at the emergence of a type of small world 
behaviour, i.e: with connecting only 10 sites the average shortest path is reduced from | to 
0.17, connecting an extra 10 sites reduces this value to 0.09. 

We see that as fat as Co < the effect of eolations is on, y to m odff y the Nations of Q 
by replacing M , the actual number of connections , with an effective one M. Expressing M 
in terms of Mq and Co alone, we find: Mq = M(l + M(q), which means that for low values 
of Co, M ~ Mq while for large values of Co the effective number of connections scales as the 
square root of the actual number of shortcuts, M ~ v/^f ■ This effect reflects the tendency 
of the shortcuts to get clustered under the influence of correlations. Hence correlations tend 
to decrease the small world effect, since the connections tend to bunch into clusters. 



III. UNDIRECTED AND CLUSTERED NETWORKS 

As far as we have N — > oo and M =finite, we can generalize our results to the cases 
where 

a) -the network has no preferred direction 
and 

b) - each site of the ring is connected to 2k of its neighbors. 

In this limit, in going from one site of the ring to another one, one travels mostly along the 
ring. Thus denoting the average shortest paths for the above cases respectively by z 3> a 
and C 2 >(, we have 

1 1 

<z> a =-<z> , <z> fe = (22) 

Z K 

from which we obtain 

Q a (z) = 2Q(2z) < z < X - (23) 
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and 

Qb(z) = kQ(kz) 0<z<-. (24) 

K 

And finally ,for the clustered undirected model one will have 

Qab(z) = 2kQ(2kz) 0<z<jj-. (25) 

IV. CONCLUSION 

We have considered the effect of correlations in a simple model of small world network, 
and shown that they generally decrease the small world effect, since under this condition the 
connections tend to bunch into clusters. More concretely in our simple model the effect of 
correlations which are controlled by a parameter (o, is to reduce (for large zetao) the actual 
number of shortcuts M to an effective one M ~ y^f > indicating a clustering of connections 
to bunches in the lattice. 

Therefore it seems that the optimal way of designing a small world network would be with 
equidistant long-range connections and in order to see the small world effect and lower the 
average shortest path, one is better to use algorithms which anticorrelate the connections. 
We have derived our results by exact analytical methods, and have shown how the transfer 
matrix technique can be used for obtaining such properties as average shortest path, or the 
distribution of shortest paths in a model of small world network. For all this we have been 
forced to study a restricted class of models. No doubt by doing computer simulations one 
can study these effects in a much broader class of models. 
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